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Source subtraction for the extended Very Small Array and 
33-GHz source count estimates 
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ABSTRACT 

We describe the source subtraction strategy and observations for the extended Very 
SmaU Array, a CMB interferometer operating at 33 GHz. A total of 453 sources 
were monitored at 33 GHz using a dedicated source subtraction basehne. 131 sources 
brighter than 20 mJy were directly subtracted from the VSA visibility data. Some 
characteristics of the subtracted sources, such as spectra and variability, are discussed. 
The 33-GHz source counts are estimated from a sample selected at 15 GHz. The 
selection of VSA fields in order to avoid bright sources introduces a bias into the 
observed counts. This bias is corrected and the resulting source count is estimated to 
be complete in the flux-density range 20-114 mJy. The 33-GHz source counts are used 
to calculate a correction to the VSA power spectrum for sources below the subtraction 
limit. 

Key words: cosmic microwave background - surveys - radio continuum: galaxies. 



1 INTRODUCTION 

Extragalactic radio sources are a major contaminant of Cos- 
mic Microwave Background (CMB) observations, particu- 
larly at small angular scales. The power spectrum {Ct = 
l{l + l)Ci/2n) of Poissonian-distributed point sources in- 
creases as £^ whereas that of the CMB decreases with in- 
creasing due to Silk damping and the incoherent addition 
of fluctuations along the line of sight through the surface of 
last scattering. The point-source foreground will therefore 
dominate CMB measurements at high £ in the absence of 
an effective source-subtraction strategy. 

Three factors complicate the implementation of schemes 
for dealing with point sources. Firstly, in the frequency 
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range ~ 30-200 GHz where many CMB experiments op- 
erate, no large-area surveys for sources exist. Simple ex- 
trapolation of flux densities from low-frequency data can 
cause worse source contamination than if no subtraction 
were attempted at all, due to the pop ulation of sources 
with rising spectra jTavlor et alJ l200lll . The positions of 
sources from low-froquoncy catalogues may be used to 
'project-out' (Bond. Jaffc & Knox 199^) sources from the 
high-frequency 30 GHz) CMB data, but at a cost of 
signal-to- noise. Secondly, 'flat '-spectrum (flux-density spec- 
tral index a ~ 0, S oc i^"*^) sources will increasingly domi- 
nate bright flux density-selected samples at higher frequen- 
cies. Flat-spectrum sources are compact objects and of- 
ten show flux density variations which means that accurate 
source subtraction at higher frequencies requires monitoring 
of source flux densities simultaneously with CMB observa- 
tions. Thirdly, in order to avoid confusion by the CMB, sur- 
veying for sources should be performed at high resolution 
where the CMB contribution is negligible. 

The Very Small Array (VSA) is a collaboration be- 
tween the Cavendish Astrophysics Group (University of 
Cambridge), the University of Manchester's Jodrell Bank 
Observatory (JBO) and the Instituto de Astrofisica de Ca- 
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narias (lAC) in Tenerife. The instrument is a 14-element in- 
terferometer operating in Ka band (26 -36 GHz) and locate d 
at Teide Observatory, Tenerife (see e.g. lWatson et aLll2003fl . 
The CMB power spectrum up to £ ~ 800 h as been mea- 
sured by the VSA in its compact configuration (Tavlor et^J 
[2003: Scott et alJl2003) and up to ^ ~ f 500 (Grainec ct^ 
I2OO3I: iDickinson et al.l f2004'l in an extended configuration. 
A unique source subtraction strategy is employed by the 
VSA to allow accurate subtraction of sources. The regions 
observed by the VSA are surveyed with the Ryle Telescope 
(RT; see e.g. lJones^ll991^ in Cambridge at 15 GHz (see e.g. 
IWaldram et al.ll2003i ~The sources detected at 15 GHz are 
then observed with a single-baseline interferometer at the 
same frequency as the VSA and at the same time. 

The purpose of this paper is to describe the source sub- 
traction observations and discuss some characteristics of the 
subtracted sources, such as spectra and variability. The 33- 
GHz source counts are estimated and used to determine 
a correction to the VSA power spectrum for unsubtracted 
sources. 



2 VSA SOURCE SUBTRACTION STRATEGY 

VSA extended array measurements currently probe angular 
multipoles up to ^ ~ 1500, where measurements of the power 
spectrum can be significantly affected by Poisson noise from 
unsubtracted radio sources. In order to mitigate this poten- 
tially serious foreground, the VSA strategy is three-fold. 

Firstly, VSA fields are selected such that there are no 
sources predicted to be brighter than 500 mjy at 30 GHz 
on the basis of the NRAO VLA 1.4-GHz All Sky Survey 
(NVS S: ICondon et al Ill 99^ and the Green Bank 4.8-GHz 
(GB6; iGreeorv et all |l996) survey. These predictions were 
made by extrapolating the fiux density of every source in the 
4.8-GHz catalogue using the spectral index between 1.4 GHz 
and 4.8 GHz. 

Secondly, the RT in Cambridge is used to survey the 
VSA fields at 15 GHz in advance of CMB observations. The 
RT synthesised beam FWHM for this purpose is 25 arcsec 
corresponding to a spherical multipole, ^ « 5 x 10'*, where 
the contribution from the CMB is negligible. The area sur- 
veyed for a typical VSA extended array field is some 13 deg'^, 
whereas the area within the FWHM of the RT primary beam 
is only 0.01 deg^. This means that a raster scanning tech- 
nique must be used so that the 15-G Hz survey can keep up 
with the VSA observations. Waldra m et al.l i2003l) describe 
the 15-GHz observations and present source lists. 

Thirdly, the sources detected at 15 GHz are monitored 
at 33 GHz simultaneously with CMB observations using a 
dedicated source subtraction baseline. 

The depth of the RT survey is determined by the 33- 
GHz source subtraction flux density limit, Sum, which in 
turn is prescribed by the requirement that the confusion 
noise from unsubtracted sources is at best very small com- 
pared with the thermal noise on a VSA map. The confusion 
noise, (Tconf, can be estimated as (■Sclieuer..l957i'l 



n{S) dS, 



(1) 



to 5'-)-d5', SI is the VSA synthesised beam solid angle and 
sources with flux densities greater than Sum have been sub- 
tracted. Then, from equation^ for a given confusion noise, 
the required Sum is given by 



•S'lim — 



f(3-/3) 



si; !^ 



1/(3-/3) 



(2) 



where the differential source count, n{S) dS, is deflned as the 
number of sources per steradian in the flux density range S 



where the differential source count is parameterised as 
n{S) = K {S/So)-". 

At the outset, the differential source counts were 
not known at the VSA observing frequency and the 
source counts from RT observations for t he Cosmic 
Anisotropy Teles c ope ( CAT: lR,obson et al.ll993h were used. 
lO'Sullivan et all (^93) found an integral count N(S > 
10 mJy) « 3 per deg^ at 15 GHz and and /3 was taken to be 
approximately 1.8. The thermal noise on a compact array 
VSA map is « 30 mJy beam~^ and the compact VSA syn- 
thesised beamwidth is ~ 30 arcmin. Assuming a conservative 
average spectral index between 15 and 34 GHz, afl — 0, the 
source subtraction limit such that the VSA compact array 
maps are thermal noise dominated is Sum ~ 109 mJy. In 
fact, a subtraction limit of 80 mJy was used for the com- 
pact array. In order that sources with rising spectral indices 
as extreme as af| = —2 are found at 5 o" at 15 GHz, the 
sensitivity of the RT survey was cr « 4 mJy. 

For the extended array observations, which measure the 
power spectrum at higher multipoles, a deeper RT survey 
is required. Using the updated counts from the compact 
array so urce monitorin g, we have n{S) = 54 S~^'^^ Jy~^ 
sr~ ( Ta vlor et al .120031) . The nominal r.m.s. noise on a VSA 
extended array map is ~ 5 mJy beam"*, the synthesised 
beamwidth is ~ 11 arcmin and the required 5iim is ~ 20 mJy 
at 33 GHz. This means that the sensitivity of the RT survey 
for the extended array must be at least ~ 0.8 mJy in order 
to find inverted-spectrum sources with af| = —2 at 5 a. 
In fact, the requirement that the RT must survey the VSA 
fields in advance of VSA observations means that a practical 
limit to the RT sensitivity is reached at 2 mJy. As a conse- 
quence, not all sources with rising spectral indices aff ^ —1 
will be found and there remains a possibility that some such 
sources exist in the VSA fields. Clearly, for measurements 
at higher £, a new strategy is required 



3 OBSERVATIONS 

The VSA source subtraction system consists of two 3.7-m 
dishes spaced 9.2 m apart, used as a single-baseline North- 
South interferometer located adjacent to the main array 
(Fig.0. The receivers, IF system and correlator are identical 
to those for the VSA. Each dish is located within an enclo- 
sure to minimize ground spill-over. With this baseline, the 
interferometer is not sensitive to large-scale CMB emission, 
while the contribution on smaller angular scales is expected 
to be negligible. Both the VSA and the source subtraction 
system measure the same combination of Stokes parameters. 
Observations of each source were made over a range of hour 
angles similar to those used by the main array. 

Over the course of the extended array monitoring pro- 
gramme, the source subtractor dishes were replaced with 
ones of higher surface accuracy and the system was upgraded 
to double sideband operation. This increased the aperture 



© 2004 RAS, MNRAS 000, [I]-?? 



Source subtraction for the extended Very Small Array and 33- GHz source count estimates 3 




Figure 1. VSA source subtracter dishes in individual enclo- Figure 2. The points indicate the locations of the 453 sources 

sures to minimise ground spill-over. The dishes are aligned North- monitored for each of the VSA mosaics. The 2° FWHM of each 

South, have diameters of 3.7 m and are separated by 9.2 m. VSA pointing is indicated by a circle. The projection is gnomonic 

and coordinates of RA{J2000) and Dec(J2000) are in degrees. 



Table 1. Nominal specifications of the VSA source subtracter. 



Location 

Maximum declination range 

No. of antennas 

Baseline separation 

Frequency range 

System temperature, Tsys (K) 

Bandwidth, Au 

Sidebands 

Correlator 

Dish diameter 

Point source flux sensitivity (l-o") 



Izana, Tenerife (2340 m) 
-5° <Dec < -1-60° 
2 

9.2 m 

26 - 36 GHz 
« 35 K 
1.5 GHz 

two 
complex 
3.7 m 
Si 0.34 Jy sl/2 



efficiency from ~ 0.5 to ~ 0.7 and doubled the effective 
bandwidth. The nominal r.m.s. of visibility amphtude for the 
source subtractor following the upgrade is cr ~ 0.34 Jy s^''^. 

The purpose of the monitoring programme is to mea- 
sure the flux densities of the sources detected by the RT at 
15 GHz at the VSA observing frequency of 33 GHz. This 
takes place simultaneously with CMB observations by the 
VSA. For a typical observation, the source subtractor drive 
system takes the coordinates of the sources within a radius 
of some 2° of the current field centre from the 15-GHz source 
list. A subset (typically ~ 30) of these sources is monitored 
over the course of each ~ 5-hour observing run, with inter- 
leaved observations of the brightest source in the field for 
phase calibration. In this manner, the flux densities of the 
sources are sampled over the course of the observations of 
a given field. The integration time for each sample is 400 
s, giving a nominal sensitivity per sample of 17 mjy and 
typically « 30 samples are obtained for each source. A to- 
tal of 453 sources were monitored for the purposes of source 
subtraction in 7 independent regions as shown in Fig. |21 



4 DATA REDUCTION AND CALIBRATION 

Data reduction for the source subtractor is performed in 
the same manner as for the main array data, using the re- 



duce software package (see e.g. IWatson et al J 120031) . The 
calibrated visibility time-series for each source in a given 
VSA field is examined by eye to check the quality of the 
phase and flux-density calibration. Bad data are flagged at 
this point and the samples for each source are averaged over 
the observing period of that fleld by the main array. 

Due to the much greater instantaneous sensitivity of 
the source subtractor compared to the VSA, the brightest 
calibrators used for the VSA (Tau A, Cas A and Jupiter) are 
not suitable for the source subtractor. In addition, Tau A 
and Cas A are resolved by the source subtraction baseline. 
The main flux density calibrator for the source subtractor is 
NGC 7027, a planetary nebula with a flux den sity of 5.45 ± 
0.2 Jy at 32 GHz, based on a measurement bv lMason et alJ 
( 1999.) . Assuming a spectral index, a = 0.1 ± 0.1, a flux 
density of 5.43 Jy at 33 GHz is adopted. The calibration 
uncertainty in the sou rce flux densi t ies is 4 per cent due to 
the uncertainty in the lMason et alj il999l) measurement. 

As with the main array, a baseline-based amplitude cali- 
bration is applied. Since the source subtractor slews between 
pointed observations, interleaved observations of phase cali- 
brator sources are made. The phase calibrator sources need 
only be bright enough to achieve good signal-to-noise ratios 
in a 400-s integration and are typically the brightest sources 
already identified in each VSA field. 



SOURCE SUBTRACTION FOR EXTENDED 
VSA DATA 



As described in lDickinson et al] (|2o3), a campaign of high- 
sensitivity measurements of the CMB up to ^ ~ 1500 was 
conducted with the VSA in an extended configuration. The 
flux densities of the sources detected at 15 GHz in the VSA 
fields were monitored by the source subtractor at 33 GHz 
during observations by the main array, as described above. 

For each field, the sources found to have mean flux 
densities over the observation period greater than 20 mJy 
were subtracted directly from the VSA visibility data. Fig.|3 
shows two examples of a VSA field before and after source 
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Table 2. The sources subtracted from the VSA6G and VSA8G 
fields in Fig. |2l The flux density measured by the source subtrac- 
tor, S, and the apparent (i.e. attenuated by the primary beam, 
as on a VSA map) flux density, Sapp, are given. 



Field 


RA{J2000) 


Dec(J2000) 


S [Jy] 


•S'app [Jy] 


VSA6G 


07^ 28™ 23" 8 


+54° 31' 16('8 


0.154 


0.115 


VSA8G 


34™ 20? 6 


+38° 57' 51^'4 


0.995 


0.089 




17^ 30™ 41! 4 


+41° 02' 56'.'6 


0.069 


0.045 




17'' 38™ IKl 


+40° 08' 19^6 


0.101 


0.056 



subtraction, demonstrating successful source subtraction 
and revealing CMB structure previously confused by the 
presence of bright sources. 

A total of 131 sources were found to have flux densities 
greater than 20 mjy and were directly subtracted from the 
visibility data. Of these 131 sources, the brightest apparent 
flux density on a VSA map was 115 mJy (see TableHJ. 

Some characteristics of the monitored sources, such as 
spectra and variability, are now discussed. 



5.1 Source spectra 

In order to examine the spectra of subtracted sources, a 
catalogue search was performed for the brightest counter- 
parts within 10 arcsec using the Astrophysical Catalogues 
Support System (CATS; iVerkho danov ot al. 1997} ■ Fig. 
shows some examples of the inverted-spectrum sources in- 
cluding J0944+3303, the only subtracted source not to have 
a positional counterpart in NVSS or FIRST, with a spectral 
index a = —2 between 15 and 33 GHz. Where available, 
data from meas urements made ove r a range of frequencies 
simultaneously iBolton et al.ll20o3 . hereafter B04) are also 
plotted. 

In order to assess whether any of the subtracted sources 
were significantly peaked, an attempt was made to fit the 
spectra to a broken power-law with the following formula 
jMarecki et al.lll999^ : 



S{u) 



^W^o)'=(l 



(3) 



where k and I are the spectral indices of the rising and falling 
parts of the spectrum respectively and So and i^q axe fitting 
parameters. In some cases, the spectrum was significantly 
peaked but a second-order polynomial of the form. 



log = So + Q log — c(log v) 



(4) 



provided a better fit to the data. Fig. |^ shows the 8 sources 
where equation |3] provided a good fit to the available data, 
or the curvature, c, of the second-order polynomial fit was 
> 0.5. Table|3|gives the associated fitting parameters. 

There is considerable interest in the identification of 
gigahertz peaked (GPS) sources and those sources with 
spectra peaking at frequencies above a few gigahertz, 
known as high frequ ency peakers (HFP). The 9C survey 
iWaldram et al.l|2003ll is the first radio-frequency survey to 
cover an appreciable area at a frequency above the 4.8- 
GHz GB6 survey with a sensitivity down to « 10 mJy. As 
such, it allows HFP candidates to be selected on the ba- 
sis of their spectra between 4.8 and 15 GHz. The follow-up 



Table 3. Spectral parameters for peaked sources of Fig. |F] The 
parameters are, peak flux density, 5maxi and corresponding fre- 
quency, fmax, the spectral indices of the rising (fc) and falling (Z) 
parts of the broken power-law fit, and the curvature, c, of the 
second-order polynomial fit. 



Name(J2000) 


'S'max [Jy] 


I'max [GHz] 


k 


I 


c 


0024+2911 


0.081 


10.2 






1.6 


0717+5231 


0.151 


2.7 






0.6 


0936+3313 


0.058 


1.8 


-0.51 


0.48 




0942+3344 


0.120 


1.1 


-0.47 


0.59 




0944+3347 


0.087 


2.3 


-0.41 


0.67 




1235+5228 


0.091 


2.5 


-0.26 


0.83 




1526+4201 


0.083 


8.3 






1.1 


1740+4348 


0.226 


7.4 


-0.48 


1.09 





observations of a subset of these sources at 33 GHz pre- 
sented here provide another high-frequency datum. An at- 
tempt has been made here to investigate if any of the sub- 
tracted sources are significantly peaked on the basis of 15- 
GHz and 33-GHz data as well as low-frequency data from 
catalogues. However, source variability means that follow-up 
observations at a range of frequencies simultaneously, such 
as those performed by B04, may be required to verify the 
existence of spectral peaks. Of the 8 sources determined here 
to be peaked-spectrum, 3 sources were also observed by B04. 
Of these 3 sources, J0024+2911 and J1526+4201 were also 
found to be peaked by B04 and their measurements of the 
third source, J0936+3313, are well fit by the falling portion 
of the broken power-law, although they lack lower-frequency 
fiux-density measurements to verify the presence of a peak 
around 1 GHz. This suggests that a high-frequency survey 
combined with a catalogue search is a useful method of se- 
lecting peaked-spectrum candidates, despite the variability 
of sources. 



5.2 Source variability 

For the majority of sources, the source subtractor lacks suffi- 
cient sensitivity for measurements of day-to-day flux-density 
variability, since we use it to measure the average flux den- 
sity over the time-scale of the VSA observations. However, 
we can get an idea of the range of variability in the sources 
monitored for source subtraction in two ways. Firstly, the 
VSA fields overlap so that the sources lying in the overlap 
are observed at different epochs. By comparing the mean 
fiux density measurements at different epochs, the variabil- 
ity of sources in the over lap regions can be examined. Sec- 
ondly, a test was applied to find those sources whose 
time-series could not be fitted by a constant fiux density. 
This latter method is only useful for bright sources observed 
frequently (e.g. phase calibrators) which thus have sufficient 
signal-to-noise. 

In order to compare the fiux densities measured at dif- 
ferent epochs, the flux densities measured at 34.1 GHz were 
interpolated to 33 GHz. Fig. 15] shows a histogram of the 
percentage variability of the 72 sources in the overlapping 
regions with fiux densities greater than 20 mJy. Of these 72 
sources, 23 were found to have varied by more than 25 per 
cent of the mean flux density, 8 by more than 50 per cent 
and 3 by more than 75 per cent. 
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Figure 3. Maps of the VSA6G (top row) and VSA8G (bottom row) fields before (left) and after (right) source subtraction. The crosses 
mark the positions of the sources predicted to have the brightest apparent flux density from source subtractor data (Table 1^. The circles 
in both maps indicate the area under the primary beam 2° FWHM. Contour levels are -30, -20, 10, 20, 30, 40, 50, 60, 70, 80, 90, 100 
and 110 mjy/beam. 



The above analysis does not take into account the er- 
ror in the flux density measurements at different epochs. In 
order to determine whether a source is variable or not, a 
test is applied, i.e. 

j=i ^ ' 

where N is the number of measurements, the Si are the in- 
dividual flux densities and AS'i their errors. This tests the 



hypothesis that the flux density time-series could be mod- 
elled by a constant. Sources for which this probability is 
0.1% are taken to be variable. Of the 453 sources mon- 
itored at 30 GHz for source subtraction purposes, only 8 
were found to be variable according to this criteriona and 
these were all phase calibrators. Fig. |7| shows the measured 
flux density time-series (as well as the 15 GHz time-series, 
where available) for the 8 sources found to be variable using 
the test. While these sources are taken to be variable for 
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Figure 4. Examples of sources with rising spectra, on tlie basis of the 15 and 33-GHz flux densities as well as the flux densities of 
positional counterparts in the CATS database (filled circles). Data from measurements by B04 are also plotted, where available (unfilled 
circles). The source J0944+3303 is the only source monitored at 33 GHz for source subtraction not to have a positional counterpart in 
NVSS or FIRST and has a spectral index of -2 between 15 and 33 GHz. 
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Figure 5. Spectra for peaked-spectrum candidate sources on the basis of the 15 and 33-GHz fiux densities as well as the flux densities of 
positional counterparts in the CATS database (filled circles). Data from measurements by B04 are also plotted, where available (unfilled 
circles). The solid curves are broken power-law or second-order polynomial fits to the data. 
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Figure 6. Histogram showing the number of sources in the over- 
lapping areas of the VSA fields whose maximum variability, ex- 
pressed as a percentage of mean flux density, falls in the ranges 
0-25, 26-50, 51-75 and 76-100 per cent. 



the purposes of further analysis, that is not to suggest that 
other sources in the sample are definitely not variable. 

For these 8 sour ces, the structure functions (see e.g. 
ISimonetti et alJ Il985il were estimated in order to check 
for the presence of characteristic timescales in the flux 
density variations. Structure function analysis allows time- 
variability to be quantified without the problems encoun- 
tered using Fourier techniques with irregularly sampled 
data. The flrst-order and second-order structure functions 
are deflned as 

D\t) = {[S{t + T)-S{t)f) 

and 

D\t) = {[S(t + 2T)-2S{t + T) + S{t)]^) 

respectively, where, S{t) is the source flux density at time, 
t, T is the 'lag' and the angle brackets denote an ensem- 
ble average. For a stationary random process, the structure 
functions are related to the variance of the process, a^, and 
its autocorrelation function, p(r), as follows: 

D\r) = 2a'[l + pir)], 

D\r) = 6a^[l + ip(2r) + ip(r)]. 

For r ^ To, where ro is the smallest correlation time-scale 
present in the series, the structure function, _D*^(r) has a 
logarithmic slope 2M. For r ^ n, where n is the great- 
est correlation time-scale present in the series, the structure 
functions reach a plateau of 2a^ for M = 1 and 6cr^ for 
M = 2. For intermediate time-scales, the logarithmic slope 
of the structure function is 5C 2M. Well-deflned structure 
functions were found for 2 of the 8 sources taken to be vari- 
able and these are plotted in Fig. |H| We now comment on 
each of these sources. 

J0958+3224 (3C232) - A compact steep-spectrum quasar 
with z — 0.53. The flrst-order structure function for this 
source has a well-defined slope of « 0.8 on time-scales of 



3-12 days, turning over at r « 12 days. The second-order 
structure function has a slope of ~ 1.3 in the same time- 
scale range, indicating that day-like structure exists in the 
time-series since purely linear trends (or segments of longer- 
term variations) would produce a flat second-order structure 
function. The time-series (Fig. shows variations of up to 
^ 30 per cent. 

J1734+3857 - This source is identified with the flat- 
spectrum BL Lac object NVSS J173420+385751 at z = 0.97. 
The well-deflned logarithmic slope of ~ 0.9 in the first-order 
structure function on time-scales of 6-39 days indicates the 
presence of structure in the time-series with a maximum cor- 
relation time-scale of ~ 39 days. The second-order structure- 
function is less clear, however the slope of « 0.9 over the 
same time-scales appears inconsistent with the presence of 
purely linear trends. The time-series (Fig. shows varia- 
tions of up to ~ 20 per cent. 

Although the test indicates variability in 8 sources, 
these are amongst the 11 brightest sources monitored for 
source subtraction and so have been observed frequently as 
phase calibrators. Fainter sources may be variable but will 
not be identified so easily as such by a test due to the 
decreased signal-to-noise ratio of these observations. Fig. |^ 
plots the ratio of observed-to-expected r.m.s. of flux density, 
for all 453 sources observed for source subtraction against 
the spectral index between 1.4 and 33 GHz, averaged in 
bins of Aqi^4 = 0.5. If differences between the observed and 
expected r.m.s. of flux density are due to variability, then 
we might expect the ratio of observed-to-expected r.m.s. to 
increase with decreasing spectral index on the basis that the 
majority of variable compact sources are flat-spectrum, and 
this appears to be the case. 



6 33-GHZ SOURCE COUNT ESTIMATES 

The set of sources subtracted from VSA data does not form 
a well-selected sample for the purposes of source-count es- 
timation. In order to estimate the 33-GHz source counts, a 
new source sample was constructed as follows. We selected 
well-defined regions within the existing VSA fields which 
had previously been surveyed at 15 GHz down to a com- 
pleteness limit of 10 mjy. The selected regions encompass a 
total area of 0.044 steradian and 370 sources with 15-GHz 
flux densities greater than 10 mJy were identified. All the 
sources identifed above the completeness limit at 15 GHz 
were observed at 33 GHz with the source subtraction base- 
line. Hence the source count we present at 33 GHz is based 
on a sample selected at 15 GHz. In addition, since the VSA 
fields were specifically chosen to have no sources brighter 
than 500 mJy at 33 GHz (on the basis of extrapolation from 
NVSS and GB6 data) we do not expect our source counts 
to be complete at high flux densities. 

We have also identifled a further bias in our data. The 
first three VSA fields (VSAl, VSA2 and VSA3, hereafter 
collectively VSA123) were chosen with an extra constraint, 
such that no source predicted to be brighter than « 250 mJy 
at 33 GHz should he within the FWHM of the primary 
beam. The later VSA fields (VSA5, VSA6, VSA7 and VSA8, 
hereafter collectively VSA5678) were chosen with the more 
relaxed point source limit of 500 mJy since their positions 
were already constrained by operational requirements. This 
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slight difference in field selection has resulted in a (deliber- 
ate) deficiency of sources with fiux densities greater than ~ 
100 mjy at 15 GHz in VSA123. This can be seen by com- 
paring the differential counts at 15 GHz for the VSA123 and 
VSA5678 fields with that determined using the 9C survey 
(see Fig. iim . In order to use sources from the VSA123 fields 
in our 33 GHz source-count estimate, we therefore attempt 
to take account of this bias. 

In using sources from VSA123 to estimate the counts 
we need to set an upper limit, 5'upper, to the flux density 
at 33 GHz to ensure that the original selection at 15 GHz 
(which is incomplete above Sis ghz ~ 100 mJy) does not 
bias the count. One approach is to examine the 15-33 GHz 
spectral indices of all sources above 100 mJy at 15 GHz in 
our whole sample. We find 12 sources, with a highest (falling) 
spectral index of 1.0. Since this a small sample, we have 
investigated the spectral index distributions for complete 
samples at lower fiux densities at 15 GHz. We find that 
approximately 80 per cent of these sources have off ^ 1.0. 
This value of 1.0 gives us an estimate for Supper of 46 mJy. 
A second approach has been to examine the 33-GHz source 
counts separately for VSA123 and VSA5678; we find the 



deficit in the VSA123 count relative to the VSA5678 count 
appears to become significant at approximately 80 mJy. We 
therefore assume that 46 mJy is a conservative estimate for 
Supper and combine the counts from VSA123 and VSA5678, 
using only those sources from VSA123 with fiux densities 
^ 46 mJy. Fig. 1111 shows the resulting differential source 
counts. The data used for the fit are shown in Table |1] A 
correction is made for the bin widths such that the number of 
sources in each bin of width, W, and centre, Sc, is multiplied 
by a factor of (1 — r^) where r = W/2Sc- 

Parameterising the differential source counts as a power 

law, 

«(S) = Ac(^) , (6) 

the least-squares fit to the power law was determined for fiux 
densities in the range 20-114 mJy. Choosing So — 70 mJy 
to lie in the centre of this range, the resulting fit gives 

/ o \ -2-34^°:- 

n(S) = (10.6l^;^) f ^^^j mJy-isr-\ (7) 

The quoted errors correspond to the 1"*^ confidence limit 
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Figure 8. The first and second-order structure functions for J0958-I-3224 (top left and riglit) and J1734-I-3857 (bottom left and right). 



Table 4. Table showing the data used to fit the differential source 
count in Fig. 1111 The normalising area for each bin is also pro- 
vided. 



Bin start 


Bin width 


Number 


Normalising 


[Jy] 


[Jy] 


in bin 


area [sr] 


0.020 


0.004 


23 


0.044 


0.024 


0.004 


20 


0.044 


0.028 


0.004 


18 


0.044 


0.032 


0.005 


14 


0.044 


0.037 


0.009 


13 


0.044 


0.046 


0.068 


14 


0.024 



for each parameter. Expressed in Jy ^ sr ^ , the source count 
is 

n(5') « 21 S--'-^^ Jy"' sr-\ (8) 

which is also plotted in Fig. 1111 The errors on the data points 
are the formal Poisson errors, in iV. 

Fig. ll2l shows the differential source count predictions of 
a 30-GHz model ( Toffolatti ct al. 1 998, hereafter T98), th e 
sourc e counts derived fr om WMAp'_jBenn«^^tal 
CBI jMason et al.ll2003l) and DASI jKovac et al.ll2002t) ex- 



periments as well as the VSA count described by equation|H| 
All counts are normalised to those expected in a Euclidean 
universe. The observed counts are consistent in the overlap- 
ping flux-density ranges, but generally fall below the T98 
model. The observed VSA counts will be used to re-scale 
the T98 model in the following section. 



7 RESIDUAL CORRECTION TO THE CMB 
POWER SPECTRUM 

The contribution to the temperature fluctuations from faint 
sources below the subtraction limit, Sum, has a component 
due to the Poissonian-distributed sources as well as a clus- 
tering term, as follows: 

rf-i2 ^ fd]3i/\ /^Poisson I ^Clusterings /^.n 

-'CMB'-'src — I Qj, I ['-^STC + '-'src ) ) 

where 

95, _ 2fc / fcTcMB \' x^e^ 

ar ^ c2 V h J (e- - 1)2 ' ^ ' 
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Figure 11. The 33-GHz differential source counts showing a fall- 
off in flux density bins < 20 mjy and > 114 mjy, ascribed to lack 
of completeness. The best-fltting power law to flux densities in 
the range 20—114 mjy (102 sources) is also plotted (solid line). 
The vertical dotted lines indicate the flux density range used in 
the flt. The vertical error bars indicate the Poisson errors in each 
bin and the horizontal error bars indicate the bin widths. 
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Figure 10. 15 GHz differential source counts derived from 
sources in VSA123 (filled circles) and VSA5678 (unfllled cir- 
cles) fields. The dashed line is the known 15 GHz source count 
(n(S) fa 51 S~^-^^) from the 9C survey. There is a (deliberate) 
deficiency of bright sources in the VSA123 fields. 



Bv is the Planck function, k is Boltzmann's constant and 
X = hv/kTcMB- We can estimate each term of equation |U] 
using (see equation 



Ci;™ = / 5' n{S) dS 
Jo 

and iScott fc Whitelll999l) 



a 



Clustering 



uJi\ I Sn{S)dS 



(11) 
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Figure 12. VSA 33-GHz differential source counts (normalised to 
those expecte d for a Euclidean un iverse) c ompared with thos e de- 
rived by CBI jMason et alJl2003^ . DASl iKovac et al.ll2002^ and 
WMAP iBennett et al.ll2003l) .The delineated regions for each 
experiment are defined by the quoted l-cr error on the normali- 
sation of n{S). The s olid curve shows the p redicted counts from 
the 30-GHz model of lToffolatti et alJ 11993) . 



where ui is the Legendre transform of the angular two-point 
function of sources, u>{9). Based on estimates of the angular 
correlation funct ion of sources from the NVSS and FIRST 
surveys (see e.g. lBlake fc Wal]||2002^ ■ we estimate that the 
contribution of the clustering term (eauation ll2ll is negligible 
over all scales measured by the VSA. 

Fig. ll3l plots the integrand of eauation llll on a log scale 
for the VSA and other experiments operating at a similar 
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Figure 1 3. S'^n( S) for the VSA compared with those derived by 
CBI ( Mason ct al. 2003), DASI iKovac e t al. 2002) and WMAP 
iBennett et al.,.2003,) . The delineated regions for each experiment 
are defined by the quoted l-tr error on the normalisation of n(S). 
T he sohd black curve sh ows the predicted counts from the model 
of lToffolatti et alJ Jl998l) . while the solid grey curve is this model 
re-scaled by a factor of 0.68 in order to best fit the VSA data. 
The source subtraction limit of 20 mjy is indicated by the ver- 
tical long-dashed line. The vertical solid line indicates the lower 
integration limit for which the residual source correction is 99 per 
cent of the total. 



frequency. The vertical long-dashed line indicates the upper 
limit of integration, 5'iim = 0.02 Jy and the vertical solid line 
indicates the lower limit of integration for which the residual 
correction is 99 per cent of the total. For differential source 
counts with power-law slopes /? < 3, the dominant contribu- 
tion to the Poisson noise comes fr om the flux decades ju st 
below the source subtraction limit jToffolatti et alJll99d) . 

The measured source counts are used to estimate the 
contribution from sources below the flux density subtrac- 
tion limit. However, it is clear that extrapolating the best- 
fit source count to lower flux densities will over-estimate the 
counts of faint sources since, as the T98 model shows, the 
counts are expected to flatten at lower flux densities. The 
VSA source counts were used to re-scale the T98 model and 
the best flt is plotted as the grey curve in Fig. ll.SI The re- 
scaling of the T98 model required to best flt the VSA 33 GHz 
source counts was found to be 0.68 ±0.07. Table|K|shows the 
residual source contribution for the VSA estimated using 
source counts measured by experiments operating at similar 
frequencies, the T98 model and that model re-scaled by 0.68 
(denoted T98*). 

The residual source correction calculated using the T98* 
model is applied to the VSA power spectrum as this method 
does not over-estimate the counts at low flux densities. 
Integrating the T98* model up to the subtraction limit 
of 20 mJy, we get T^mbCsic = 130 x lO"'^ sr (or, 

207 at ^ = 1000, see Table EJ. This means that the 

residual source power spectrum is estimated to be AT^rc = 
207 X (^/lOOO)^ ^.K^ and this is binned and subtracted di- 
rectly from the CMB band-power estimates as shown in 
Fig.m 



The T98 30 GHz model must be re-scale d by 0.75 in or- 
der t o fit the Ka-band WMAP source counts iArgiieso et alJ 
l2003h . and this is consistent at the 1-a level with the re- 
scaling of 0.68 ± 0.07 required by the VSA source counts. 

The cosmological parameter estimation incorp orating 
the extended-array power spectrum dat a of^Dickinson et all 
(2004!) was performed bv iRebolo et alJ i2004ft . In order to 
asses s the effectiven ess of our source subtraction procedure, 
Rcb olo et alj (|200^ included the parameter Ax in the like- 
lihood analysis, where Ce = 2tvAx x 10~®. A significant 
AT^ _ Ax£^ component of the power spectrum may be 
considered as an indication of an error in our estimate of 
the residual source correction assuming that the contribu- 
tion from the SZ effect is negligible for £ ^ 1500. 

The combined WMAP and VSA datasets were used 
for this analysis as well as two external priors, the Hubble 
Space Telescope Key Project constraint o n Hubble's con- 
stant of Ho = 72 ± 8 km sec~^ Mpc~^ iFreedman et alJ 
200 ll and constraints on lar ge-scale structur e from the 2dF 
Galaxy Redshift Survey JColless et alJl200J) . The analysis 
gives Ax = -46 ± 132 and Ax = -86 ± 123 with the HST 
and 2dF priors respectively. In both cases, A^ is consistent 
with zero and assuming that point sources are the only sig- 
nificant AT^ (X component, the source subtraction has 
been successful to at least within ~ 100 ]iK^. 



8 IMPLICATIONS FOR OTHER CMB 
EXPERIMENTS 

iMason et alJ (l2003l) have previously found power in excess 
of the expected primary anisotropy in the multipole range 
^=2000-3500 of the TT power spectrum from an analysis 
of deep field observations with the CBI in 2000. The results 
of further TT measure ments have recently been announced 
jReadhead et al.ll2004l) which, when combined with previ- 
ous data, act to reduce the signiflcance of the excess power 
detection. In this section we examine the implications of the 
33-GHz source monitoring programme for the original CBI 
dee p-field detecti o n of ex cess power. 

iMason et alj i2003l) report a band-averaged power 
508 ^K2 in their highest £ bin (2010 < £ < 3500) which 
is 3.1-u in excess of their band-averaged model power of 
66 \xK^. The detection of excess power could be due to a 
population of sources with strongly inverted spectra not ap- 
pearing in NVSS and therefore not 'projected-out' of the 
CBI data or appearing in the residual source correction. All 
sources detected at 15 GHz and subtracted from the VSA 
fields had positional counterparts in NVSS except one, which 
had a strongly inverted spectrum {a = —2) between 15 and 
33 GHz. There was no evidence, therefore, of a new popu- 
lation of strongly inverted sources at 15 GHz. However it is 
important to note that, with a lower completeness limit of 
20 mJy at 33 GHz, the source count estimate is insensitive 
to sources with rising spectral indices less than a « — 1. A 
deeper survey will be required to investigate whether there 
is a significant population of sources appearing at 33 GHz 
with spectra more inverted than this. 

The excess power could also result from an underesti- 
mate of the residual source correction due to sources be- 
low the NVSS detection threshold of 3.4 mJy. This cor- 
rection is estimated to be C""" = 0.08 ± 0.04 Jy^ sr~^ by 
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Table 5. Residual source power for a subtraction limit of 20 mjy calculated using the source counts (n(5) = k Jy~^ sr"^) measured 
by the VSA and other experiments as well as the T98 30-GHz model and that model re-scaled by 0.68 (T98*). The expected residual 
source contribution at £ = 1000 is also given. Due to the expected flattening of the counts at lower flux densities, source counts measured 
at high flux densities (e.g. WMAP) will over-estimate the residual source power for subtraction to low flux densities. 
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Figure 14. The final, double-binned, VSA extended array power spectrum before (unfilled circles) and after (filled circles) the estimated 
residual source power spectrum (dotted line) has been subtracted. 



iMason et alJ l)200,^) using Monte-Carlo simulations of source 
populations based on the 31 GHz source counts determined 
from CBI deep and mosaic maps {n{S) = 92 Jy~^ sr"^) 
and the OVRO-NVSS distribution of spectral indices with 
— 0.45. In order to account for the expected flattening 
of the counts, the re-scaled T98 model was used to deter- 
mine the residual contribution with the 1.4 GHz limiting flux 
density of 3.4 m Jy extrapolated to 33 GHz using the mean 
NVSS-VSA spectral index af^^^ = 0.54 observed by the VSA. 
This mean spectral index was derived from the 31 sources 
in VSA5678 fields with 533 ghz ^ 30 mJy, in order to ensure 
that only high signal-to-noise 33-GHz measurements were 
included. This results in an estimate of 6"^"° — 0.03 Jy^ sr~ ^ 
which is consistent with that found bv lMason et alJ i2003l) . 
It is unlikely, therefore, that Poissonian-distributed point 
sources constitute the dominant contribution to the CBI 
excess. However, it is important to note that the observed 
NVSS-VSA spectral index distribution is subject to various 
selection effects (discussed in Section |SJ which will be ad- 
dressed by source-subtraction schemes for higher-resolution 
VSA observations. 



9 SOURCE SUBTRACTION FOR THE 
SUPER-EXTENDED VSA 

The VSA has previously observed in configurations probing 
£ up to 1500. In principle, observations up to higher £ are 
possible, limited by the size of the tilt-table upon which the 
horn are mounted. The array has now been reconfigured to 
probe £ up to 2500. In order to maintain the filling factor of 
this 'super-extended' array, the 322 mm mirrors have been 
replaced with lightweight 550 mm mirrors. The front-end 
amplifiers have also been upgraded, resulting in a system 
temperature decrease from 35 to 28 K. 

Based on the T98* model, the required source subtrac- 
tion limit such that the residual source power spectrum is 
less than a model CMB power spectrum at ^ = 2500 is 
Slim ~ 2 mJy at 33 GHz, as shown in Fig. 1151 This deeper 
source subtraction at 33 GHz also has implications for the 
15-GHz survey with the RT. In order that all sources with 
rising spectra a = —2 are found at 5-a, the required comple- 
tion limit at 15 GHz is around 0.4 mJy, requiring a RT sen- 
sitivity of ~ 80 ti.Jy. Planned upgrades to the RT as part of 
the Arcminute Microwave Imager programme iKneissl et alJ 
I2OOL) will result in an increase in fiux sensitivity by a factor 
of ~ 10 and hence in survey speed by ~ 100, making surveys 
of the required depth feasible. 
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Figure 15. The residual source power spectrum for a source sub- 
traction limit of 10 mjy (dashed-dot line), 5 mjy (dashed) and 
2 mjy (dotted line). Also shown is a model CMB power spectrum. 
A source subtraction level of ~ 2 mjy is required for a resid- 
ual source power spectrum less than the CMB power spectrum 
at ^ = 2500. 



An alternative approach involves the use of an instru- 
ment currently under development . The goal of the One Cen- 
timetre Receiver Array (OCRA; iBrowne et alJ l2000l) pro- 
gramme is to develop a 100-beam receiving system for all- 
sky surveys at 30 GHz. The first stage (OCRA-p) involves 
building a 2-beam prototype and this has been installed on 
the Torun 32 m dish in Poland for testing. The first astro- 
nomical task of OCRA-p will be to make an unbiased survey 
of the discrete sources in the VSA fields down to a few mJy. 
The brightest sources detected using OCRA-p could then 
be monitored for variability by the VSA source subtracter 
during CMB observations by the main array. 



10 CONCLUSIONS 

The VSA source subtracter has been used to monitor at 
33 GHz sources in the VSA fields detected by the Ryle Tele- 
scope at 15 GHz. A total of 131 sources with fiux densities 
greater than 20 mJy were subtracted from the data. 

A test indicated that 8 of the subtracted sources were 
variable and two of these had well-defined structure func- 
tions showing evidence of variability on timescales of days. 
The source subtracter does net have sufficient instantaneous 
sensitivity to confirm variability in the fainter sources, how- 
ever there is evidence of a correlation between the ratio of 
observed-to-expected r.m.s. in the time-ordered flux density 
data and the spectral index between 1.4 and 33 GHz. The 
subtracted sources include five new peaked-spectrum candi- 
date sources, with estimated peak frequencies in the range 
1-7 GHz. 

Based on areas within the VSA fields surveyed by the 
Ryle Telescope to 10-mJy completeness, the 33-GHz differ- 
ential source counts were estimated to be n{S) ~ 21 
Jy~^ sr~^. That is, there are estimated to be 0.9 sources per 
square degree brighter than 20 mJy. The 33-GHz counts pre- 



sented here are estimated to be complete in the ffux density 
range 20-114 mJy. 

Previous VSA analyses have extrapolated the source 
count to fiux densities below the source subtraction limit 
in order to estimate the residual source contribution. This 
method takes no account of the expected ffattening of the 
counts at faint fiux densities and so is likely to over-estimate 
the residual contribution. Instea d, the 30 G Hz differential 
source- count model of [joffolatti e t al.l Jl99d) was re-scaled 
by 0.68 to best fit the VSA data. Integrating the re-scaled 
T98 model up to 20 mJy provided an estimate of the resid- 
ual source contribution and this was subtracted from the 
estimated power spectrum. The residual contribution was 
estimated to be ^ 207 H-K^ at ^ = 1000. The re-scaling 
of the T98 model required by the VSA and WMAP source 
counts indi cate that the as s umpt ions about source spectra 
adopted bv lToffolatti et all l|l99^ may need to be refined. 

Using the re-scaled T98 model and the 1.4-33 GHz spec- 
tral index statistics, th e residual source co rrection for the 
CBI deep-field data of iMason etTaO J2003l) was esimated. 
We find a contribution of 0.03 3y^ sr~^ which is consistent 
with t he value of 0.08 ±0.04 Jy^ sr"^ found bv lMason et alJ 
J2003h . 
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APPENDIX A: 33-GHZ SOURCE LIST 

In this section, the 33-GHz source list used for the source 
count estimates is presented. For each source found to have 
a ffux density greater than 20 mJy at 33 GHz, the measured 
ffux density, S, and 1-a thermal error, AS, in Jy are pre- 
sented in Table Al. In addition to the quoted error, there is 
a calibration uncertainty of ~ 4 per cent in fiux density. 
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Table Al. 33-GHz source flux densities. 



Name(J2000) RA(J2000) Doc(J2000) S [mjy] AS [mjy] 



0010+3044 


00'' 


10"" 


08?4 


30° 


44' 


42" 


20.0 


3.0 


0010+2717 


00'' 


10"" 


28! 7 


27° 


17' 


54" 


25.8 


6.1 


0011+2803 


00'' 


11"" 


33! 8 


28° 


03' 


46" 


29.4 


6.0 


0011+2928 


00'' 


11"" 


46! 


29° 


28' 


28" 


27.3 


3.0 


0012+2702 


00'' 


12"' 


38! 1 


27° 


02' 


40" 


27.5 


6.3 


0012+3053 


00'' 


12"" 


50! 3 


30° 


53' 


19" 


31.5 


3.3 


0014+2815 


00'' 


14m 


33! 8 


28° 


15' 


08" 


34.8 


3.7 


0018+2907 


00'' 


18"' 


50!9 


29° 


07' 


40" 


22.9 


2.3 


0019+2817 


00'' 


19"' 


08! 8 


28° 


17' 


56" 


23.9 


2.4 


0019+2956 


00'' 


19"' 


37!8 


29° 


56' 


02" 


30.9 


2.3 


0023+2928 


00'' 


23™ 


39! 7 


29° 


28' 


22" 


35.5 


2.6 


0024+2911 


00'' 


24™ 


34! 5 


29° 


11' 


29" 


29.8 


2.5 


0024+2724 


00'' 


24m 


50! 5 


27° 


24' 


27" 


29.7 


4.0 


0028+3103 


00'' 


28™ 


11!2 


31° 


03' 


35" 


23.0 


5.0 


0028+2914 


00'' 


28™ 


17!1 


29° 


14' 


29" 


30.2 


4.0 


0028+2954 


00'' 


28™ 


31!9 


29° 


54' 


52" 


22.4 


3.2 


0029+2958 


00'' 


29™ 


40!3 


29° 


58' 


51" 


22.0 


3.0 


0259+2434 


02'' 


59™ 


10!5 


24° 


34' 


12" 


20.0 


5.1 


0301+2504 


03'' 


01™ 


43! 9 


25° 


04' 


43" 


25.3 


3.4 


0301+2442 


03'' 


01™ 


46! 1 


24° 


42' 


58" 


25.0 


5.2 


0303+2531 


03'' 


03™ 


59! 5 


25° 


31' 


35" 


30.0 


2.9 


0307+2656 


03'' 


07™ 


07!7 


26° 


56' 


07" 


20.6 


2.9 


0309+2738 


03'' 


09™ 


22! 


27° 


38' 


56" 


113.0 


3.2 


0310+2810 


03'' 


10™ 


18!8 


28° 


10' 


02" 


28.9 


5.1 


0311+2530 


03'' 


11™ 


40!9 


25° 


30' 


12" 


55.2 


5.0 


0313+2718 


03'' 


13™ 


44! 9 


27° 


18' 


48" 


30.0 


6.1 


0314+2636 


03'' 


14™ 


58! 3 


26° 


36' 


49" 


41.0 


4.0 


0708+5409 


07'' 


08™ 


30!5 


54° 


09' 


05" 


20.7 


8.4 


0712+5430 


07'' 


12™ 


44! 6 


54° 


30' 


27" 


23.8 


5.9 


0714+5343 


07'' 


14™ 


47!4 


53° 


43' 


27" 


29.8 


5.9 


0715+5609 


07'' 


15™ 


16!8 


56° 


09' 


55" 


22.9 


8.5 


0716+5323 


07'' 


16™ 


40!4 


53° 


23' 


13" 


110.6 


5.9 


0717+5231 


07'' 


17™ 


30!0 


52° 


31' 


04" 


31.3 


8.9 


0719+5229 


07'' 


19™ 


16!4 


52° 


29' 


56" 


26.6 


8.8 


0719+5526 


07'' 


19™ 


37!9 


55° 


26' 


07" 


34.6 


6.2 


0720+5213 


07'' 


20™ 


38! 8 


52° 


13' 


30" 


29.2 


8.8 


0723+5355 


07'' 


23™ 


22! 6 


53° 


55' 


38" 


35.6 


3.3 


0724+5653 


07'' 


24™ 


20! 7 


56° 


53' 


43" 


20.9 


13.3 


0727+5621 


07'' 


27™ 


45!4 


56° 


21' 


14" 


30.2 


12.9 


0728+5325 


07'' 


28™ 


03! 3 


53° 


25' 


15" 


52.9 


3.8 


0728+5431 


07'' 


28™ 


23! 8 


54° 


31' 


17" 


112.4 


3.8 


0731+5338 


07'' 


31™ 


16!9 


53° 


38' 


57" 


39.9 


3.6 


0733+5605 


07'' 


33™ 


28! 7 


56° 


05' 


42" 


51.5 


12.5 


0739+5343 


07'' 


39™ 


43! 6 


53° 


43' 


07" 


23.2 


4.8 


0931+3049 


09'' 


31™ 


03! 1 


30° 


49' 


12" 


32.6 


9.5 


0934+3050 


09'' 


34™ 


47! 3 


30° 


50' 


52" 


21.9 


3.6 


0936+3313 


09'' 


36™ 


09! 5 


33° 


13' 


09" 


23.8 


3.4 


0942+3344 


09'' 


42™ 


36! 2 


33° 


44' 


37" 


25.0 


4.0 


0944+3115 


09'' 


44™ 


11!6 


31° 


15' 


24" 


25.3 


3.6 


0944+3347 


09'' 


44™ 


20! 1 


33° 


47' 


56" 


34.8 


3.8 


0946+3050 


09'' 


46™ 


13!5 


30° 


50' 


22" 


28.5 


3.9 


0950+3201 


09'' 


50™ 


48! 7 


32° 


01' 


43" 


21.1 


5.0 


1215+5349 


12'' 


15™ 


05! 9 


53° 


49' 


55" 


46.0 


5.0 


1215+5336 


12'' 


15™ 


28! 8 


53° 


36' 


08" 


114.0 


16.0 


1215+5154 


12'' 


15™ 


46! 1 


51° 


54' 


50" 


26.0 


6.0 


1216+5244 


12'' 


16™ 


23! 6 


52° 


44' 


01" 


20.0 


10.0 


1219+5408 


12** 


19™ 


43! 3 


54° 


08' 


33" 


29.0 


10.0 


1221+5429 


12'' 


21™ 


18!2 


54° 


29' 


06" 


20.7 


3.5 


1223+5409 


12'' 


23™ 


13! 1 


54° 


09' 


09" 


34.9 


3.5 


1229+5522 


12'' 


29™ 


09! 3 


55° 


22' 


30" 


109.2 


8.4 


1229+5147 


12'' 


29™ 


22! 


51° 


47' 


07" 


41.3 


4.9 


1234+5054 


12'' 


34™ 


16!0 


50° 


54' 


23" 


33.2 


5.4 


1235+5228 


12'' 


35™ 


30! 6 


52° 


28' 


28" 


27.9 


2.6 


1235+5340 


12'' 


35™ 


48! 


53° 


40' 


05" 


62.1 


3.5 
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Table Al - continued 33-GHz source flux densities. 



Name(J2000) RA(J2000) Dec(J2000) S [mjy] AS [mjy] 



1237+5057 


121^ 


37m 


23! 6 


50° 


57' 


21" 


79.0 


6.0 


1238+5325 


12*1 


38" 


07!6 


53° 


25' 


54" 


27.2 


3.6 


1240+5334 


12*1 


40™ 


05? 5 


53° 


34' 


33" 


27.1 


6.2 


1240+5441 


12^ 


40™ 


14? 2 


54° 


41' 


48" 


27.8 


8.8 


1240+5334 


12^ 


40™ 


42?4 


53° 


34' 


23" 


24.8 


7.0 


1241+5141 


12^ 


41™ 


16?5 


51° 


41' 


29" 


22.0 


5.0 


1523+4156 


15^ 


23™ 


09?4 


41° 


56' 


25" 


27.9 


4.3 


1523+4420 


15^ 


23™ 


34?4 


44° 


20' 


31" 


23.9 


10.9 


1526+4201 


15^ 


26™ 


45? 3 


42° 


01' 


41" 


34.2 


3.9 


1528+4219 


15*1 


28™ 


00? 1 


42° 


19' 


14" 


27.4 


2.5 


1533+4107 


15*^ 


33™ 


27?9 


41° 


07' 


23" 


28.1 


3.3 


1539+4217 


15*1 


39™ 


25? 7 


42° 


17' 


27" 


34.0 


2.3 


1539+4123 


15'' 


39™ 


36? 8 


41° 


23' 


33" 


33.7 


3.3 


1541+4114 


151^ 


41™ 


01?2 


41° 


14' 


28" 


32.5 


3.2 


1545+4130 


15*1 


45™ 


21?4 


41° 


30' 


27" 


30.2 


4.0 


1546+4449 


15^ 


46™ 


04? 6 


44° 


49' 


13" 


24.0 


2.0 


1546+4257 


15^ 


46™ 


22? 5 


42° 


57' 


56" 


23.6 


4.1 


1723+4206 


Ijh 


23™ 


25? 1 


42° 


06' 


10" 


24.4 


3.4 


1726+3957 




26™ 


32? 6 


39° 


57' 


02" 


109.5 


6.3 


1727+4221 


Ijh 


27™ 


49? 3 


42° 


21' 


41" 


44.3 


2.6 


1728+3935 


17^ 


28™ 


32? 1 


39° 


35' 


16" 


24.1 


6.3 


1730+4102 


17^ 


30™ 


41?4 


41° 


02' 


57" 


39.3 


2.7 


1731+4101 


17^ 


31™ 


23? 6 


41° 


01' 


38" 


25.8 


2.2 


1733+4034 


lYh 


33™ 


07?5 


40° 


34' 


28" 


20.8 


6.2 


1734+3857 




34™ 


20?5 


38° 


57' 


52" 


939.3 


1.9 


1737+3958 




37™ 


41? 1 


39° 


58' 


18" 


36.1 


6.4 


1737+4154 


lYh 


37™ 


59?4 


41° 


54' 


54" 


28.1 


2.4 


1738+4008 


17^ 


38™ 


19? 1 


40° 


08' 


20" 


105.7 


6.4 


1738+4026 


17^ 


38™ 


49? 7 


40° 


26' 


21" 


20.6 


6.4 


1740+4234 


Ijh 


40™ 


51?9 


42° 


34' 


47" 


57.0 


5.0 


1744+4014 




44™ 


25? 1 


40° 


14' 


49" 


97.3 


6.1 


1745+4059 




45™ 


28? 5 


40° 


59' 


52" 


267.7 


6.2 
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